On the formation of Wigner molecules in small quantum dots 
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It was recently argued that in small quantum dots the electrons could crystallize at much higher 
densities than in the infinite two-dimensional electron gas. We compare predictions that the onset 
of spin polarization and the formation of Wigner molecules occurs at a density parameter ~ 4 
to the results of a straight-forward diagonalization of the Hamiltonian matrix. 
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I. INTRODUCTION 



If the number of electrons artificially confined on a 
quasi two-dimensional electron island (made for example 
in a semiconductor heterostructure) is very large, many 
properties of such a so-called "quantum dot" or "arti- 
ficial atom" [| can be described from what is known 
about the limit of the infinite (two-dimensional) electron 
gas (2DEG). Until now, most experiments were per- 
formed at electron densities which are slightly below the 
equilibrium density of the 2DEG. The liquid-like proper- 
ties then still dominate. For only a few trapped particles 
(as experimentally realized in vertical quantum dots [^), 
pronounced addition energy maxima as a consequence of 
shell structure and aligned spins in the mid-shell regions 
due to Hund's rules were observed in close analogy to 
atomic physics ||^. Even the simplest picture of N non- 
interacting particles in a two-dimensional harmonic trap 
could explain many features of the conductance spectra. 
For larger systems, mean field approaches like Hartree- 
Fock [IP^ or density functional methods ||6|-[l^ were ap- 
plied. In the small- limit, much theoretical work fo- 
cused on exact diagonalization techniques ]r^ , [T^ . This 
approach was mostly used for dots in magnetic fields, 
where correlations become increasingly important with 
field strengths. It was particularly successful in the (in- 
teger and fractional) quantum Hall regime where one can 
restrict the basis set to the lowest spin-polarized Lan- 
dau level Hjie). Quantum Monte Carlo (QMC) ll7|-|2^] 
methods provide alternative approaches yielding energies 
whose accuracy reaches that of exact diagonalizations. 

When the electron density is lowered and the Coulomb 
energy increases relative to the kinetic energy, correla- 
tions begin to strongly dominate the electronic struc- 
ture also in the absence of magnetic fields. For densi- 
ties smaller than a certain critical value a Wigner crys- 
tal 1^ will be formed, in which the Coulomb inter- 
action distributes the single electrons classically on a 
lattice. For the homogeneous two-dimensional electron 
gas, such crystallization is expected at very low densi- 
ties. Monte-Carlo calculations indicate that in the 2D 
bulk a transition to a Wigner crystal-like state, pre- 
ceeded by a transition to a polarized phase [g2|, oc- 
curs only at densities corresponding to Wigner-Seitz radii 
rs,2D > 37 ajj j2^ (the density np and rs,2D are re- 
lated by no — l/(7''^s 2Li))' whereas in 3D the classical 
limit lies as high as Vs^sd = 100 [|^. (In the fol- 



lowing, for simplicity we write rs 2D — fs-) Chui and 
Tanatar |^ found that in 2D systems without trans- 
lational invariancc the critical density parameter for a 
fiuid-solid transition is shifted to a considerably smaller 
value {vs ~ 7.5 a*g). Could this be important for finite 
systems such as the above mentioned lateral or verti- 
cal semiconductor quantum dot structures? This ques- 
tion was recently posed [p| p^j2^ and it was argued that 
in finite systems confining only a few particles localiza- 
tion would indeed occur at significantly higher densities 
than in the 2D bulk. In the Wigner limit the few elec- 
trons in the trap would distribute such that their elec- 
trostatic repulsion is minimized. The internal structure 
of the wave function of the many body system should 
then have the symmetry of the corresponding classical 
charge distribution. Wigner crystallization was found to 
be particularly pronounced in quantum dots with steep 
walls and polygonal geometry jS^I . Egger et al. have 
performed quantum Monte Carlo studies using a mul- 
tilevel blocking algorithm |^^. For parabolic quantum 
dots with azimuthal symmetry they reported that at a 
critical density of =4 a*g the formation of Wigner 
molecule-like ground states should become energetically 
favorable. Hartree-Fock (HE) methods by definition fail 
to accurately describe the correlated regime. When per- 
formed in an unrestricted scheme, however, spontaneous 
symmetry breaking and localization in the spatial distri- 
bution of the electronic densities of quantum dots and 
lateral quantum dot molecules at rg 
tributed to the onset of Wigner crystallization |^ . 

In the present article we report numerically exact con- 
figuration interaction calculations. This method has a 
long history in quantum chemistry, and was applied to 
quantum dots by many different groups ]l^ . Much of 
the previous work, however, concentrated on the elec- 
tronic structure in large magnetic fields where the elec- 
tron gas is polarized. Our purpose here is a comparison of 
the exact diagonalization results to the above mentioned 
recent predictions of localized states in the low-density 
limit and zero magnetic field. We first give a brief outline 
of the configuration interaction method and then turn to 
a discussion of the many-body spectra of a six-electron 
quantum dot at zero angular momentum as a function 
of the average electron density in the dot. Calculations 
for different r^-values indicate that the ground state re- 
mains unpolarized. At values of which are accessible to 
exact diagonalization techniques, for a dot confining six 
electrons clear signals of formation of a Wigner molecule 
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could not be observed. We conclude with a brief com- 
parison to results of density functional theory (DFT) in 
the local density approximation (LSDA) describing the 
electronic ground state structures. 



II. METHOD AND CONVERGENCE 

Consider N interacting electrons trapped in a circu- 
larly symmetric harmonic well V{r) — m* Lo^r^ , where 
= -|- 2/^. (In the quasi two-dimensional limit one 
assumes that the confinement in z-direction is much 
stronger than in the x-y-plane. Then, only the lowest 
subband in z-direction is populated.) We write for the 
Hamiltonian 
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Here, to* and e are the effective mass and the dielectric 
constant. The calculations are done for different values 
of the density parameter rs which determines the aver- 
age particle density in the dot, uq = l/(7rr^). The lat- 
ter is approximated by setting the oscillator parameter 
uJq ~ e^/(47reoeTOrf a/TV) 0- Throughout this paper we 
use effective atomic units in which the length unit a*g is 
a factor e/m* times the Bohr radius qb, and the energy 
is given in effective Hartree, Ha* — Ha-(TO*/£^). (For 
GaAs for example, m* = 0.067TOe and e = 12 A, for which 
the length and energy units then scale to a*g — 97.9 A 
and Ha* = 11.9 meV.) To diagonalize the Hamiltonian, 
Eq. (|l|), the spatial single-particle states of the Pock space 
are chosen to be eigenstates of the two-dimensional har- 
monic oscillator with optimized oscillator parameter uj. 
In general, the electron-electron interaction tends to ex- 
pand the system, and thus lu ^ loq. This effect becomes 
stronger with increasing and we used the (empirical) 
relation uj — loq/^/tI . The two-body matrix elements 
of the electron-electron interaction are calculated using 
the addition theorem for 1/ | ri — r2 | To set up 

the Fock states for diagonalization, we use eight lowest 
oscillator shells containing 36 states and sample over the 
full space with a fixed number of spin down and spin up 
electrons, Ni — N. From this sampling, only those 

Fock states with a given total orbital angular momentum 
and a configuration energy (corresponding to the sum of 
occupied single-particle energies) less than or equal to a 
specified cutoff energy Ec are included (see Figure |l| be- 
low for an example) . The purpose was to select only the 
most important Fock states from the full basis, hereby 
reducing the matrix dimension to a size d ~ 10^. To ob- 
tain all the eigenstates we have to set — N-^ — N/2 
for even particle numbers {Sz — 0, all states with dif- 
ferent total spin have this component), and analogously 
we would have Ni 



iV| ± 1 for odd numbers. Once 
the active Fock states are at disposal, the Hamiltonian 
matrix is calculated. For Lanczos diagonalization we use 
the Arpack library [E9|. Finally, the total spin of each 



eigenvector is determined by calculating the expectation 
value of the operator. 

As mentioned above, setting up the Hamiltonian ma- 
trix from chosen Fock states and subsequent diagonal- 
ization only in principle yields an exact solution of the 
many-body problem. For reasons of numerical feasibility 
it is necessary to truncate the set of basis functions to be 
used in the diagonalization. One then has to make sure 
that convergence of the spectra is reached with respect to 
the cut-off. As the required matrix size increases rapidly 
with iV, computational expenses severely restrict the cal- 
culations to only the smallest systems at not too large 
values of r^. Thus, with increasing electron number or 
Vs , the results become less accurate due to the restricted 
number of basis states that can be included in the cal- 
culations. (The fact that the ground and excited states 
are the closer in energy the larger the particle number N 
imposes an additional difficulty.) 
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FIG. 1. Convergence of the many-body spectra of 
iV = 6 electrons in a harmonic trap at a density corre- 
sponding to = 4 a*g as a function of the cut-off energy 
Ec/{Tllu). Shown are the 6 lowest states with L — 0, 
together with the spin S of each state. 



For a quantum dot confining = 6 electrons at a den- 
sity corresponding to = 4 (the largest value of rs 
we found accessible within the calculational scheme used 
here) , Figure |^ shows the convergence of the many-body 
spectra as a function of the cut-off energy Ec- The low- 
est possible Fock state for six electrons has two particles 
in the state \n,l) = |0,0) and four particles in |0, ±1). 
Thus, the configuration energy equals 2huj + 4:-2huj. This 
means that for the spectra with different cut-off ener- 
gies displayed in Figure |l| all excitations up to an energy 
Ec — lOhuj are included. The many-body spectra for lAfko 
and I6huj differ drastically from the results obtained for 
Ec > 20?iu!. Looking at the relative ordering of the levels 
and the spin sequence, it becomes clear that convergence 
is only reached for Ec > 20fiuj. The ground-state energy 
for zero angular momentum is Eq — 3.049 Ha* for 22Tilu 
and Eq — 3.045 Ha* for 2AhLU. An extrapolation to in- 
finite cut-off energy can be made by plotting the total 
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energy as a function of {E^ — 10?iaj)~'^/^ |3^. This gives 
the estimate 3.043 Ha* for the fully converged results at 
Ts = 4 a|j. For Ec — lAhw a too small number of Slater 
determinants was included to build up the required cor- 
relations, such that the polarized S — 3 state appeared 
as the ground state. (We identify a similar effect in unre- 
stricted HF results mentioned above Q , where the single 
Slater determinant that is available incorrectly favors a 
spin-polarized ground state Q.) While for Ec — 22huj 
the matrix dimension 44181 with 21448811 non-zero ma- 
trix elements is reasonably small, the value E^ = 2ATilo 
already yields a matrix dimension 108375 with 67521121 
non-zero matrix elements. As for larger densities the 
states are less correlated, a smaller number of Slater de- 
terminants is needed for an accurate description. Den- 
sity parameters larger than = 4a^ or a higher number 
of particles than iV = 6 would go beyond the limits of 
numerical feasibility and accurate results could not be 
obtained. 



III. MANY-BODY SPECTRA OF A 
SIX-ELECTRON QUANTUM DOT 

We now analyze the many-body spectra and the se- 
quence of spins for the low-lying states as the two- 
dimensional density parameter is varied. We choose 
the particle number TV = 6 as it corresponds to the small- 
est dot size for which classically two stable crystalline 
structures co-exist: A pentagonal ring with one electron 
at the center, and a slightly distorted six-fold ring pl| , p2| . 
We fix the angular momentum to L = and show in Fig- 
ure ^ the 50 lowest states for a quantum dot confining 
N = 6 particles at density parameters between rs — I a*^ 
and = 4 a|j (in steps of 0.5 a*^). To obtain a better res- 
olution of the spectra the energies of the eigenstates Ci are 
scaled such that the energy difference between the ground 
state and 50*'^ excited state equals one, i.e. plotted are 
the dimensionless quantities q = (cj — ei)/(e5o — ei) . 
(The total energies of the ground state with = and 
the excited state with = 3 are given in Table 1 below.) 
At a very large density corresponding to 1 a|j, the 
ground state has spin S* = and is separated from the 
lowest excited state with spin 5 = 2 by a gap of 0.49 Ha*. 
This state is followed by a state with S ~ 1 and again 
another (excited) spin singlet. The lowest fully polarized 
state is only found at a fairly high energy, the energy dif- 
ference to the S = ground state being about 0.95 Ha*. 
(We note that for A'' = 6 at = 1.73 a*^ we obtained ex- 
cellent agreement of the S = ground state energy with 
the result of Pederiva et al. |]l^.) As rg is increased, the 
fully polarized S* = 3 state moves down in energy. At 
Vg = 2.5 a*g it has passed the excited singlet, but is still 
far from competing with the nonpolarized S = ground 
state. From the evolution of the energy difference of the 
S = 3 state to the ground state, we do not expect any 
crossing of the polarized state and the ground state un- 
less Ts becomes much larger than 4 a*^ . Estimating the 
decrease in energy of the polarized state with respect to 



the ground state as is increased, our data seem to sup- 
port the result of Egger et al. that the ground state 
of the 6 electron dot is not polarized for Tg-values smaller 
than about 8 a*n. 



FIG. 2. Spectra for angular momentum L = and six 
electrons in a harmonic trap at densities corresponding 
to Wigner-Seitz radii between Vg = I a*g and = 4 ajj 
in steps of 0.5 a*g. The square and the triangle show the 
first excited singlet state and the lowest fully polarized 
state, respectively. The lowest state has always S — 0. 
The energy axis is scaled such that the energy difference 
between the ground state ei and 50*^* excited state £50 
equals one, i.e. plotted are the dimensionless quantities 
Ei = (fi - ei)/(e50 - ei) 

Table 1 compares the total energies of the ground state 
with spin zero and the lowest polarized state with spin 
5 = 3 with the corresponding result obtained from DFT, 
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TABLE 1: Table of energies [Ha*] for the paramagnetic 
{S — 0) and ferromagnetic {S = 3) states in a 6-electron 
quantum dot for different densities. In the diagonaliza- 
tion we used Ec = 22fiLu for < 3.5 a*g and Ec = 2AhLo 
for Ts = 4 a|j. For comparison the energies obtained 
with the local spin density approximation (LSDA) are 
also shown. 




3 







Tg — oQj Q 






T — An^ 

1 g '±U, Q 






plppfrnn inf prj^r'l'irnic; frpj^fpH in T jS1T~) A l-<^nr fViP T^T-'^'X' 

^I^-i-CL- Li VJJ.i J-llLCJ. CXV^ LlWllO ID Li CCXLC-VJ. Ill U kJ J— ' -Ti . i \JL ullC- l-Zi 1 


E [Ha*] 


Q 






Q 




rpmilts wp impH an intpmolatinn fovrniil?! fov tlip TanfitpiT'- 

1 Lll UO VV LIO v3 vl CA'll 111 U v^l vJlCt Ul vJll 1\_J1 111 LllClr 1\_J1 Ull^I^ 1 Ci'llCi' UCll 










r^p-pip-plp^;" ^9*^1 PYpVifinCTP-pnrrplfii'inn pnprcrv AA/p rf^^f^r NTil 

v^CL/v^in^y tJl t-^^jiidiit'^ ijVJi J. t,ici uivjii ciici ti y . vvc i cici liyj I'l 




n 
u 


n 
u 




n 
u 


n 
u 


Tny i"np iiTrl'nPT' ripi"ailQ fon^Pi-mno' "hnp nnmpripfll mp1"nnrl 

HJl LilH^ ILAlUllC-l vH^IjOjIID VjVjllV.-Vjl 11111 tl LllO IIU.111'^1 IL-Clil lllt-LlUJU.. 




1 

j_ 


1 


O.UtJ4: 


1 


1 


Tlip r^T pnprp'ips of flip S — PTonnrl ^tafp^ and S — ^ i^n- 

11"^ \_yJ- •^llV'l ^I'-'O KJl. UllV^ k-J KJ ^1 LillvJ. OLCIu^-jO Cdllvl 1^ t J lOv-' 


A 1 Q4 


1 
± 


■? 

o 


OfiO 
o .uuu 


9 


n 


iTiPv rTinrnflvp tstpII Tsrifn fnp T,m^A tpqiiI^'Q T^nv "t"np r^flvam— 

1111^1 l^UlilUcXi C Well WILjII LliO J_JO-L/iT- iCoLiiljO. i Ui LllO Uclri OjIII 


A 1 


9 


n 

u 


O.UuZ 


1 
1 


Q 
O 


afnptir' rasp tViP TiSTlA civps Inwpr pnprp'ips than tViP px- 

Clf^lie L Ivj VjCA^O j U llVj 1 J 1 y 1 \. ^1 V VjO Iv-* VV Vjl V^llV-vl ^IV-yO UllCAll UllV' V-'jV 


4 901 


1 
± 


1 
± 


o .uuo 


1 
± 


1 

± 


cict results, when Vg ^ 2 ci^- This might be ni8,inly due to 


4 90^ 







O .UU?J 





1 


thp fart that thp Tanatar-r^pnprlpv intprnnlatinn formnla 

Ullv^ IClVj U UllCIf U Ullvj X clllcl L CAjI v-' L'v-'l Ivj y lilU^l l_/ V^lcl L IWll IWl 111 Lllcl 


A 900 


1 
± 


n 

u 


o .uuu 


1 
± 


n 


slightly overcstinicitGS the correlation energy. Nevertlie- 


4 90Q 



u 


9 


o .uuo 


9 


1 

± 


less, the LSDA gives surprisingly accurately the energy 


4.209 





3 


3.068 





2 


difference between the fully polarized (5" = 3) and the 


4.213 


1 


2 


3.070 


1 


2 


paramagnetic (5 = 0) state, as seen in Fig. (For com- 


4.216 


2 


1 


3.070 





3 


parison, we also show the results for the infinite electron 


4.216 


2 


2 


3.071 


3 





gas.) 


4.225 


3 












0,3' 



0) 

cd 
E 

CD 
k_ 

to 

Q. 
LU 



0,2' 



■Je 0,1 



o 

Q. 



0,0' 



infinite gas 
6-dot exact 
6-dot LSDA 



0,0 



0,6 



1/rs 



FIG. 3. Energy between the fully polarized (5 — 3) 
and paramagnetic state {S = 0) for six electrons as a 
function of l/r^. The result for the infinite 2D gas (per 
6 electrons) is calculated using the interpolation formula 
of Tanatar and Ceperley |E3| . 



Figure ^ shows that for = 4 a|j there are only two 
L = states between the ferromagnetic {S = 3) state 
and the paramagnetic (5 — 0) ground state. However, if 
we consider all L- values, there are several states within 
this energy range. This is shown in Table 2 where ener- 
gies, spins and angular momenta of all levels up to the 
first ferromagnetic state are given. Indeed, the lowest 
excited state has L = 1 and S = 1. It is possible that at 
very large Tg this partially polarized state might become 
the ground state instead of the fully polarized S = i 
state in. 



TABLE 2: Table of energies [Ha*], spins S and angular 
momenta L of all levels up to the lowest ferromagnetic 
state for = 3a^ {left) and Ts = Aa*g [right). 



IV. CHARGE DENSITIES AND PAIR 
CORRELATION 



For = 4 ajj the radial densities of the 5 = ground 
state and the 3rd excited state of L = 0, which is the 
lowest fully polarized state with spin 5 = 3, are shown in 
Figure |4| The azimuthally symmetric charge density for 
the polarized case shows a clear maximum at the center 
surrounded by an outer ring of lower density. In the para- 
magnetic case, the density profile is more smooth, the 
maximum density being at about r « 6 a|j. The LSDA 
result shows a clear minimum at the origin, while the ex- 
act result has a larger density at the center. For compar- 
ison with the results of Egger et al. [|l^ we also show the 
density p{r) multiplied by a factor 27rr (cf. lower panel 
of Fig. ^). For the polarized state, the maximum at the 
center is now seen as a clear shoulder in the density pro- 
file. Note that this is missing in the paramagnetic ground 
state density. (This is in disagreement with the results of 
Egger et al. [|l8| who found for the ground state a density 
profile with a clear shoulder as in the polarized case.) The 
azimuthal averages of the density profiles qualitatively 
have similarities with the broken symmetry solutions of 
the unrestricted HF [^ which for the paramagnetic case 
results in a ring of six electrons, and for the ferromagnetic 
case (ground state in HF) a ring of five electrons with one 
electron in the center. However, the localization of the 
electrons is largely exaggerated in HF. Opposite to HF, 
the LSDA correctly gives the paramagnetic state as the 
ground state, and its density profile resembles the exact 
result. LSDA does not break the azimuthal symmetry 
until > 8a*g when spin- or charge density wave-like 
states can occur [§ . Purely classical Monte Carlo [||j3ll 
computations have shown that for < 6 the charges 
are distributed on the perimeter of the dot, and none of 
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the particles occupies the dot center. This changes for 
N = 6, where the charge distribution with lowest energy 
consists of five electrons sitting on a ring, with the re- 
maining electron occupying the center of the dot. This 
configuration is labeled by (5,1). If all 6 particles are 
arranged on the dot perimeter (labeled by (6,0)), the 
classical state is stable but has a higher energy than the 
(5, l)-configuration. 
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FIG. 4. Charge density of a dot confining = 6 elec- 
trons in a harmonic trap at = 4 ajj; the exact result 
{solid line) is compared to the LSDA result {dashed line). 
Shown is the density n{r) (upper panel) and n(r) multi- 
plied by 27rr (lower panel) for the paramagnetic 5* = 
state (left, H) and the ferromagnetic 5 = 3 state (right, 
'11). The S = 3 state is separated from the S = ground 
state by 0.026 Ha*. 



The classical charge distribution can be arbitrarily ori- 
ented. The density from the CI solution, however, must 
be circularly symmetric. For an azimuthal average of 
the (5, l)-pentagon structure, one would expect a pro- 
nounced maximum of the electron density in the center, 
and a less pronounced maximum at the dot radius. Cor- 
respondingly, the (6, 0) configuration should correspond 
to a minimum of charge density in the center and a max- 
imum at finite radius. A first comparison of exact di- 
agonalization calculations with the results of the mean- 
field approximation was given by Pfannkuche et al. p^ ] 
for "quantum dot helium", i.e., quantum dots contain- 
ing only two electrons. They found from a compari- 
son of exact diagonalizations with Hartree- and Hartree- 
Fock results that the exchange and correlation contri- 
butions are crucial. While the triplet state showed a 
reasonable agreement between the exact result and HF, 
the singlet could not be well reproduced. As mentioned 
above, Yannouleas and Landman |^ reported that in ge- 
ometrically unrestricted HF at a density corresponding 
to « 3.5 a|j the A'^ = 6 ground state is polarized and 
shows enhanced localization in the charge density. This 
5 = 3 state exhibits the same geometry than the clas- 



sical distribution of 6 electrons in a harmonic well: 5 
particles are equidistanly localized on the perimeter of 
the dot, and the 6th particle is trapped in the center of 
the harmonic well. The non-polarized 5 = state cor- 
responding to the (6, 0)-configuration is about 0.034 Ha* 
higher in energy. The exact diagonalization results de- 
scribed above do not support these HF results. Although 
being limited to small rg values due to the necessary re- 
strictions of the basis set, the systematic evolution and 
energy sequence of CI energies and densities shown in 
Figs. 1^ and ^ seems to indicate that polarization as well 
as formation of Wigner molecules in circularly symmet- 
ric, parabolic wells would be impossible at densities as 
large as predicted by HF |^. The geometrically unre- 
stricted solution of the Kohn-Sham equations of tend 
to overestimate the r^-value at which spontaneously bro- 
ken spin- or charge-symmetries can occur in the internal 
structure of the wave function |33|. Although calculated 
in a geometrically unrestricted DFT scheme, the fully 
converged LSDA densities for rs < 4a|j shown in Fig. ^ 
are azimuthally symmetric. Although LSDA suffers from 
the self-interaction problem, at the densities in question 
the results are in better agreement with CI studies than 
the unrestricted HF results. 
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FIG. 5. Pair correlation functions g-\a{4>) calculated at 
the outer maxima of the density distribution (cf. upper 
panel in Fig. ^) as a function of the angle 4> between the 
electrons for the ground state {lower panel) and the ex- 
cited polarized state with 5 = 3 {upper panel). The insets 
show schematic pictures of the electron configurations. 



It is finally of interest to study how the spin and spatial 
symmetry in the internal structure of the wave function 
can be recognized in the pair correlation function 



9]<y{'^) = {n^{r,Q)n„{r,ip)) 



(2) 
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which describes the probabhty to find another (spin up 
or down) particle if a particle with spin up is placed at 
(r, 0) . Here r is the radius of maximum density and Lp 
the angle between the electrons. Figure || shows gi;i{^) 
and (?||('^) for the ground state (lower panel) and g^^{(p) 
for the excited polarized state with 5 = 3 (upper panel). 
From the 93- values of the maxima in g^a{'^) one clearly 
concludes that the 5 = state has 6-fold symmetry with 
antiferromagetic spin ordering, whereas the fully polar- 
ized case shows four maxima, corresponding to a five-fold 
symmetry. These intrinsic symmetries are in qualitative 
agreement with the unrestricted HF results, although the 
crystallization predicted by HF does not yet occur. 



V. CONCLUSIONS 



We commented on the recent conjecture that Wigner 
molecules would form in quantum dots at rather large 
electron densities ^ 3.5 Pp8|. Our results are 
essentially exact up to < Aa^ for six confined parti- 
cles. The many-body spectra, densities and pair correla- 
tions obtained for = 6 clearly illustrate that the on- 
set of formation of Wigner molecules and in particular, 
polarization of the ground state should be expected at 
much higher r^-values than anticipated from unrestricted 
HF 1^ . The critical density at which such a transition oc- 
curs does, in fact, strongly depend also on geometry ||2^ ] 
and on the number of confined particles at fixed aver- 
age electron density in the dot. For the six-electron dot 
at Ts ~ 4a|j in question the "exact" ground-state 
clearly prefers 5 = and shows antiferromagnetic order 
in the pair correlation. The polarized state with C5„- 
symmetry is clearly higher in energy even than the 5 = 1 
state, which would also be a candidate for the classical 
(5, 1) ground state configuration |3^,^l[]. In addition, for 
values below = A a*g we did not find clear signals 
of rotational structure in the spectra for non-zero angu- 
lar momenta that would indicate a crystallized ground 
state. We note that the situation is different for N < 5, 
where for densities as large as = 2ag the low-lying 
states could be well understood by assuming a square- 
shaped (4, 0) Wigner molecule for the internal structure 
of the wave function and analyzing its rotational struc- 
ture. This also became clear when comparing the low- 
energy spectrum of a Heisenberg model with four elec- 
trons on a square [Q. For N > 5, however, this simple 
picture does seem to not hold, as our results for = 6 
clearly point out. The ground-state energies and densi- 
ties obtained by density functional calculations in the lo- 
cal spin density approximation agree rather well with the 
results of exact diagonalizations, even though this com- 
parison was restricted to a small particle number where 
the accuracy of the local density approximation is ques- 
tionable. This gives some confidence that the method 
is well suited for describing the ground state electronic 
structures for larger sizes. 
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